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FOREWORD 


This  work  is  an  interim  report  on  studies  of  a  fundamental  nature  concerning 
the  transfer  of  energy  from  shock  waves  into  the  internal  modes  of  molecules  in 
solid  explosives  and  has  application  to  the  study  of  initiation  and  shock 
sensitivity  of  explosives.  The  work  is  supported  principally  by  Air  Force 
Office  of  Scientific  Research  Project/Task  2301/A6,  AFOSR-MIPR-82-00004. 
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CHAPTER  1 
INTRODUCTION 


OBJECTIVE 

This  is  an  interim  report  on  work  which,  it  is  hoped,  will  lead  to  a 
greater  understanding  of  how  the  energy  from  a  shock  wave  is  transferred  to  the 
internal  vibrational  modes  of  molecules  in  solid  explosives.  Thus,  the  work 
should  be  applicable  to  the  study  of  the  initiation  and  shock  sensitivity  of 
explosives.  If  we  could,  for  example,  predict  shock  induced  energy  transfer 
rates  for  particular  types  of  bonds,  then  we  could  design  energetic  materials 
with  preferred  initiation  and  sensitivity  characteristics. 


HISTORY 

In  1979,  D.  J.  Pastine  and  co-workers  pointed  out  that  because  many 
explosives  have  weak  intermolecular  bonds  and  strong  covalent  intramolecular 
bonds  the  average  frequency  of  intramolecular  vibrations  (internal  mode  or 
optical  mode  vibrations),  (u0,  would  be  much  higher  than  the  average  frequency  of 
intermolecular  (lattice  mode  or  acoustic  mode)  vibrations,  ioa,  with  the  ratio 
u0/iua  typically  of  the  order  of  ten.^  Thus,  the  immediate  effect  of  a  shock  on 
such  materials  would  be  to  increase  the  temperature  of  the  acoustic  vibrational 
branch  while  leaving  the  optical  branches  at  the  initial  temperature.  The 
relaxation  time  required  before  the  internal  molecular  temperature  reaches  a 
critical  value  sufficient  for  the  shock  to  grow  to  detonation  would  be 
sufficiently  long  so  as  to  be  comparable  to  and  indeed  determine  the  shock  pulse 
duration  required  to  produce  detonation  at  a  given  shock  pressure.  In  other 
words,  the  relaxation  time  for  thermal  equilibration  of  the  internal  modes  is 
the  controlling  factor  in  the  initiation  of  reactions.  Using  a  simple  classical 
mass  and  spring  model,  Pastine  et  al.  estimated  that,  at  a  shock  temperature  of 
500°K,  the  lower  limit  to  the  acoustic/optical  relaxation  time  is  of  the  order  of 
several  microseconds  for  systems  in  which  the  acoustic  frequencies  (ua  are  around 
1013  rad/sec  and  the  optical  mode  frequencies  iu0  are  around  10^  rad/sec.  They 


^-Pastine,  D.  J.  ,  Edwards,  D.  J.  ,  Jones,  H.  D.  ,  Richmond,  C.  T.  ,  and  Kim,  K.  , 
"Some  New  Concepts  Relating  to  the  Initiation  and  Failure  of  Detonable 
Explosives,"  in  High  Pressure  Science  and  Technology,  Vol.  2,  ed  by 
K.  D.  Timmerhaus  and  M.  S.  Barber,  Plenum,  New  York,  1979. 
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also  predicted  that  the  relaxation  times  would  be  very  sensitive  functions  of 
the  frequency  ratio  w0/wa  with  the  relaxation  times  decreasing  by  four  orders 
of  magnitude  when  the  ratio  decreases  from  10  to  6.  Thus,  the  very  lowest 
frequency  intramolecular  modes,  which  are  typically  bending  rather  than 
stretching  modes,  should  be  the  most  important  in  the  relaxation  process. 

The  work  of  Pastine,  et  al. ,  was  followed  by  that  of  E.  T.  Toton  who 
developed  a  more  refined  quantum  mechanical  description  which  exploited  the 
disparity  between  intra-  and  inter-molecular  mode  vibrational  frequencies.*  The 
results  of  Toton' s  calculations  also  pointed  toward  the  importance  of  the  lowest 
frequency  intramolecular  modes  in  the  relaxation  process.  In  this  work,  we  will 
develop  Toton's  model,  calculate  shock  induced  internal  mode  transition  rates 
for  nitromethane ,  and  attempt  to  relate  our  results  to  critical  shock  initiation 
data  for  nitromethane. 


CHARACTERIZATION  OF  THE  SHOCK 

When  a  shock  wave  travels  through  a  solid  material,  it  excites  the  lattice 
normal  vibrational  modes  to  higher  levels  or,  in  other  words,  creates  acoustic 
phonons.  The  shock  also  compresses  the  solid,  increasing  the  frequency  of  the 
acoustic  modes.  We  wish  to  calculate  the  rate  at  which  acoustic  mode  energy  is 
transferred  to  the  molecular  internal  modes  (relaxation  rate).  For  this 
purpose,  it  is  assumed  that  the  distribution  of  acoustic  mode  energy  relaxes  to 
a  thermal  distribution  in  a  time  which  is  short  compared  to  the  time  required  to 
create  a  significant  number  of  optical  (internal  mode)  phonons.  Thus,  one 
immediate  effect  of  the  shock  is  to  raise  the  temperature  of  the  acoustic  modes 
while  leaving  the  internal  modes  "cold."  The  non-radiative  transition  rates 
between  internal  mode  levels  then  give  an  estimate  of  the  rate  at  which  the 
internal  modes  relax  to  the  new,  higher  temperature.  Later,  we  will  need  to 
determine  the  change  in  the  average  acoustic  mode  frequency  produced  by  the 
shock.  This  can  be  deduced  approximately  from  the  compression  by  integrating 
the  Gruneisen  parameter  along  the  shock  Hugoniot.  Hence,  in  this  model,  the 
effect  of  the  shock  can  be  characterized  by  two  quantities,  the  compression  and 
the  acoustic  mode  temperature  produced  by  the  shock. 


ACOUSTIC  MODES  THERMALIZE  ON  A  PICOSECOND  TIME  SCALE 

The  assumption  of  rapid  thermalization  of  the  acoustic  modes  compared  to 
the  internal  modes  is  not  unreasonable  since  internal  mode  frequencies  are 
generally  much  higher  than  lattice  mode  frequencies  and  this  mismatch  should 
greatly  reduce  the  rate  of  energy  transfer  between  lattice  and  internal  modes. 
From  the  work  of  Van  Vleck,  we  can  estimate  that  the  acoustic  modes  should 
thermalize  on  a  picosecond  time  scale. ^  Van  Vleck  computed  the  rate  of  energy 
transfer  between  lattice  oscillators  due  to  anharmonic  perturbations  when 
different  portions  of  the  frequency  spectrum  are  not  in  thermal  equilibrium. 


*Toton,  E.  T. ,  NSWC/WO,  unpublished. 

^Van  Vleck,  J.  H. ,  "Calculation  of  Energy  Exchange  between  Lattice 
Oscillators,"  Physical  Review,  59 ,  730  (1941). 
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His  result  for  the  transition  rate  from  level  n  to  level  n-1  for  a  particular 
lattice  oscillator  when  the  temperature  of  the  main  body  of  lattice  oscillations 
is  at  a  temperature  T  is 

„  5 

Vn-1  'v'  nD(T/TD)2  (l) 

where  vD  is  the  Debye  frequency,  is  the  Debye  temperature,  and  D  is  a 
»  constant  characteristic  of  the  solid.  The  value  of  D  is  about  10“^  sec^ 

for  typical  solids.  Equation  (1)  is  valid  for  T  >  Tq.  For  n=l ,  T=300°K, 
and  vq  ~  3x10^2  Hz,  Equation  (1)  gives  transition  lifetimes  of  the  order 
of  10--*-2  seconds. 


BORN-OPPENHEIMER  APPROXIMATION  APPLICABLE  TO  INTERNAL  MODES 


Van  Vleck  used  the  usual  first  order  perturbation  theory  to  calculate  the 
transition  rate.  While  this  is  satisfactory  in  the  case  of  lattice  relaxation, 
it  will  prove  unsatisfactory,  in  most  cases,  for  lattice-internal  mode  relaxation 
since  internal  mode  frequencies  may  be  an  order  of  magnitude  greater  than 
lattice  mode  frequencies.  An  exchange  of  energy  between  lattice  and  internal 
modes  which  creates,  for  example,  one  optical  phonon  will  involve  the 
annihilation  of  five,  ten,  or  more  acoustic  phonons.  To  determine  the 
transition  rate  in  this  case  would  require  a  fifth,  tenth,  or  greater  order 
perturbation  calculation.  Thus,  Toton*  developed  a  Born-Oppenheimer  formalism 
analogous  to  the  adiabatic  approximation  used  in  treating  the  coupling  of 
localized  electronic  states  to  lattice  vibrations^  and  similar  to  a  formalism 
described  by  Lin.^  In  this  formalism,  the  lattice  modes  are  considered  the 
slow  subsystem  components  and  the  internal  modes  are  considered  the  fast 
subsystem  components.  The  approximation  gets  better  as  the  disparity  between 
internal  and  lattice  frequencies  increases. 


NON-AD I ABATIC  OPERATOR  INDUCES  NON-RADIATIVE  INTERNAL  MODE  TRANSITIONS 

In  the  Born-Oppenheimer  approximation,  the  total  Hamiltonian  H  of  the 
system  is  separated  into  two  parts,  the  adiabatic  part  K  and  the  non-adiabatic 
part <£.  The  Born-Oppenheimer  states  are  eigen-solutions  of  K.  The  non-adiabatic 
part  can  be  considered  an  interaction  which  induces  non-radiative  transitions 
between  the  stationary  Born-Oppenheimer  states. 


vToton,  E.  T. ,  NSWC,  unpublished. 

^Perlin,  Yu.  E. ,  "Modern  Methods  in  the  Theory  of  Many-Phonon  Processes," 
Soviet  Physics-Uspekhi,  6,  542  (1963). 

^Lin,  S.  H.,  "Theory  of  Vibrational  Relaxation  and  Infrared  Absorption  in 
Condensed  Media,"  J.  Chem.  Phys.,  65,  1053  (1976). 
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AN  APPROXIMATE  EXPRESSION  FOR  THE  NON-RADI ATIVE  TRANSITION  RATE 

In  order  to  obtain  an  expresion  for  the  non-radiat ive  transition  rate 
between  internal  mode  levels,  we  expand  the  interatomic  potential  in  a  series  of 
terms  up  to  the  third  order  in  the  vibrational  normal  coordinates.  The  order  of 
magnitude  of  the  third  order  coefficients  is  obtained  from  dimensional  arguments 
to  be  described  later.  In  addition,  we  make  an  approximation  of  the  Condon  type 
in  which  we  assume  that  certain  internal  mode  matrix  elements  do  not  depend  on 
acoustic  mode  coordinates.  We  also  assume  that  the  acoustic  modes  can  be 
characterized  by  a  single  average  frequency  which  we  will  refer  to  as  the  Debye 
frequency.  Then,  isolating  a  single  internal  mode  and  ignoring  the  interaction 
between  internal  modes,  we  obtain  (from  the  non-adiabatic  operator)  the 
following  rate  for  a  transition  from  internal  mode  level  n  to  level  n+1. 

(2n+l)[2n+l  +  (2n+5)to  /to  ]+4m^/a)^ 

<Wl  -  <n*l>  - £_jd - W01  (2) 

1  +  5u)n/(i)  +4u)2/uj2 
o  Do 


W01  =  4-^n*  (1  +  5a)£)/w0 


2  2  /u>  . 

+  4o)  /to 

D  O  "(Dn  \T* 


(00/U)] 


-fio)0/2kT 

e 


(2a) 


X* 


fiu)  2(0,, 
o  _ D 

ke  Aft 


(2b) 


In  Equation  (2a),  T  is  the  temperature  of  the  thermalized  acoustic  modes,  u>0 
is  the  circular  frequency  of  the  internal  mode,  and  ioD  is  the  (circular) 

Debye  frequency.  The  anharmonic  coupling  effects  are  all  condensed  into  the 
quantity  An,  called  the  Stokes  shift.  Equations  (2),  (2a),  (2b)  are  valid  for 
temperatures  T  greater  than  the  Debye  temperature  TD  (=fia)p/k,  k  is  Boltzmann's 
constant)  but  less  than  T$oD/An.  We  will  see  later  than  <oD/An  »  1.  For  the 
Born-Oppenheimer  approximation  to  be  valid,  it  is  required  that  Wq/idj  >>  1 .  To 
obtain  the  transition  rate  Wn  change  the  sign  of  the  argument  of  the 

exponential  function  in  Equation  (2a)  and  change  n+1  to  n  in  the  first  factor  in 
Equation  (2). 


CHARACTERISTICS  OF  THE  NON- RADI AT IVE  TRANSITION  RATE 

Several  inferences  may  be  made  from  Equation  (2)  and  (2a).  For  transitions 
between  low  lying  levels,  transitions  between  the  ground  and  first  excited 
states  will  have  the  longest  lifetimes.  The  internal  modes  with  the  lowest 
frequencies  u  generally  will  have  the  greatest  non-radiative  transition 
rates  due  to  ?he  factor 

“0/“D 

(T/T*) 

This  factor  also  can  lead  to  large  increases  in  transition  rate  for  relatively 
small  increases  in  temperature  and  Debye  frequency. 
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INTERPRETATION  OF  THE  STOKES  PARAMETER 

The  Stokes  parameter  Afi  which  appears  in  Equation  (2a)  has  a  simple 
interpretation  for  optically  active  internal  modes.  The  coupling  between 
internal  and  lattice  modes  splits  a  single  internal  mode  absorption  or  emission 
line  into  a  band  of  lines  and  the  maximum  for  the  emission  band  is  lower  in 
frequency  than  the  maximum  for  the  absorption  band  by  an  amount  Afl .  In 
addition,  the  second  moment  of  the  band  is  a  function  of  Afi.^  If  the  band 
shape  were  gaussian,  there  then  would  bt  a  straightforward  relation  to  the  band 
width.  This  is  the  case  when  the  formalism  is  applied  to  electronic  states 
(F-centers  in  alkali  halides).  However,  as  we  will  see,  the  band  shape  for 
internal  vibrational  modes  is  highly  non-gaussian  so  that  it  is  difficult  to 
extract  A 8  from  the  band  width. 

In  the  next  chapter  we  discuss  the  formal  development  of  the  theory. 


->Lax,  M.  J.,  J.  Chem.  Phys.,  20,  1752  (1952). 
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CHAPTER  2 


THEORY 


USE  BORN-OPPENHEIMER  APPROXIMATION 

We  apply  the  Born-Oppenheimer  approximation  to  the  normal  vibrational  modes 
of  a  molecular  lattice.  The  internal  modes  constitute  the  fast  subsystem  and 
the  acoustic  modes  constitute  the  slow  subsystem.  The  separation  of  the  total 
Hamiltonian  into  an  adiabatic  part  and  a  non-adiabatic  part  X  which  induces 
non-radiative  transitions  between  Born-Oppenheimer  states  is  described  in 
Appendix  A. 


GOLDEN  RULE  TRANSITION  RATE  EXPRESSION 

Let  n^  be  the  acoustic  mode  quantum  number  for  a  mode  of  frequency  .  Let 
i  be  an  internal  mode  quantum  number  for  an  initial  internal  mode  state  and  let 
f  be  the  quantum  number  for  the  final  internal  mode  state.  Let  n  denote  the  set 
of  initial  acoustic  mode  quantum  numbers  {%}  and  let  m  denote  the  set  of  final 
acoustic  mode  quantum  numbers  (m^j.  Then  the  transition  rate  from  the  initial 
to  final  state  produced  by  the  non-adiabatic  interaction  is  obtained  from  first 
order  time-dependent  perturbation  theory  (golden  rule). 

W(i,n-f,m)  =  2it  |  <f>ra  |«C|i,n>  |  2  6(Efm  -  Ein)  (3) 

T1 

We  are  interested  in  the  total  transition  rate  assuming  a  thermalized 
distribution  of  acoustic  mode  levels,  so  we  must  sum  Equation  (3)  over  final 
acoustic  states  m  and  average  it  over  initial  acoustic  states  n.  The  result  is 


W(i-f)  =|l  £  Pn|<f,m|<£|i,n>l2  6  (Efm  -  Ein) 


where  Pn  is  the  probability  that  an  acoustic  state  with  quantum  numbers  {n^}  is 
realized.  For  a  thermalized  distribution,  this  is  a  Boltzmann  probability 
distribution. 


Pn  =  Q_ie  n 


Z-(3E 
e  » 
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6  =  1/kT 


E 


n 


'y  *  ^n<+^ 


(5c ) 
(5d ) 


The  validity  of  using  the  golden  rule,  Equation  (3),  for  the  transition  rate  has 
been  studied  by  Lin  using  a  master  equation  approach. ^ 


CONDON  APPROXIMATION 

The  remainder  of  this  discussion  will  be  concerned  with  evaluating  the 
transition  rate  given  by  Equation  (4).  When  the  slow  subsystem  is  a  collection 
of  harmonic  oscillators  the  form  of  the  non-adiabatic  operator  £  is  given  in 
Appendix  B,  Equation  (B-2).  Using  this  expression,  we  obtain  the  matrix  element 

,  .  j,,  v—'  (  .  a<t>  ■  ,  a ^ . 

<fm|<£|  in>  =  -  y.  ^ic-kMfml— i  — — >  +  ±  <*f*fm|  — —■ 

<  (  39k  3q<  1 

If  we  assume  that  the  reduced  matrix  elements 


.  34<  • 

dCcfi  ■  -  .He 

2 

yr\Ka  -  - 1  **>,c 

2  3  q| 


(7a) 


(7b) 


do  not  depend  on  the  acoustic  coordinates  q*  (our  "Condon  approximation"), 
then  we  may  write 


<^fm  |«£|  in)> 


E  k«i- 

tr  ' 


3  fyy.  | 

fm|- — ^— ■>  +  wKfi  <4ifmkin>) 
3  9ic  ) 


(8) 


While  the  Condon  approximation  is  often  a  good  one  in  dealing  with  electronic 
states,  its  validity  is  more  questionable  in  this  case.  However,  the 
simplification  it  produces  is  considerable.  Similarly,  it  is  often  assumed  that 
the^-fi  term  in  Equation  (8)  can  be  neglected  compared  to  the^j^  term  since 
theT^f^  term  contains  the  second  derivative  of  the  wave  function  with 
respect  to  the  slow  system  coordinates.  Because  of  the  questionable  validity  of 
this  assumption  and  because  dropping  the  term  leads  to  no  appreciable 
simplification,  we  will  retain  it. 


^Lin,  S.  H. ,  "On  the  Master  Equation  Approach  of  Vibrational  Relaxation  in 
Condensed  Media,"  J .  Chem.  Phys . ,  61  ,  3810  (1974). 
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INTEGRAL  EXPRESSION  FOR  THE  TRANSITION  RATE 


The  sum  appearing  in  Equation  (4)  can  be  converted  to  an  integral  by 
several  methods  summarized  by  Perlin  in  his  review  article  (See  Reference  3, 
footnote  p.  3).  One  convenient  method  involving  the  use  of  the  integral 
representation  of  the  delta  function  is  outlined  in  Appendix  C.  The  resulting 
transition  rate  (see  Equation  (C-9))  is 


W(i+f ) 


=  i_  J  Ffi(t)exp|iuifit  +T:  F<(t)A|cfi  -  s|dt 


(9) 


where 


Ffi(t)  =  k  f  i  I^k  ( c )  +  ^  'j£<  f  iAx  f  i^ic  ^ c  ^  ^  1  <£*X  f  iAx  f  i%  ^ 

K  K  < 

+  ^  ^  ^icfi^Xfi  +  fi^X  fi^AX  f  i^X  ^ 


K  ,X 


(9a) 


1 


cosh(iu)  t  +  Mi™  ) 
F^<t)  -  k  2 


2  sinh  (i  Sfuiv ) 


(9b) 


Elf(t)  “  FK(t)  -  i  coth  (i  Sfuo^ ) 


(9c) 


S  -  ^  ^  i  coth  (i  Bfia)(C  )  4Kfi 


(9d) 


Ak  f  i  f  “  9Ki 


(9e ) 


In  Equation  (9e),  qKf  and  qK£  are  the  equilibrium  values  of  the  dimensionless 
acoustic  normal  coordinate  q  when  the  internal  mode  is  in  states  f  and  i, 
respectively.  We  define  our  dimensionless  coordinates  to  be  those  which  put 
the  kinetic  energy  operator  in  the  form  shown  in  Equation  (B-l)  in  Appendix  B. 


See  Appendix  F. 


•1 

ft 
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APPROXIMATE  EVALUATION  OF  THE  TRANSITION  RATE  INTEGRAL 

An  approximate  evaluation  of  the  integral  in  Equation  (9)  may  be  obtained 
by  the  method  of  steepest  descents  as  described  by  Perlin  (Reference  3, 
footnote,  p.  3)  and  outlined  in  Appendix  D.  The  result  is  still  somewhat 
unweildy  so  we  make  the  further  approximation  that  the  sums  over  the  acoustic 
modes  can  be  replaced  by  sums  in  which  the  frequencies  are  replaced  by 
some  average  frequency  For  convenience,  we  will  refer  to  as  the 

Debye  frequency  although  the  connection  is  somewhat  tenuous.  However,  there  is 
some  evidence  that  the  average  frequency  determined  from  the  lattice  infrared 
absorption  spectrum  correlates  well  with  the  Debye  frequency  determined  by  other 
methods,  for  example,  from  specific  heat  data.^ 

The  transition  rate  expression  resulting  from  these  approximations 
reproduced  below. 

W(i-*-f)  =  i_  JIE —  Ffi(-iz0)  exp{g(z0)}  (10) 

*2  Kg"(z0) 


where 


g"(z0)  =  t/p 


exp{g(z0)}  =  expjl  fifia)Dp+P-S^ 


IPI 


Ffi(-iz0)  -  'y  ^  |  f  j  |  +  |  fiV  fi|2|FD-S/2S0(2 

IC  IC  V  ) 

*  2  XXfAfij  |fD  -  S/2S0| 

♦EX  fii2 

K 


S  *  SQ  coth 


P  =-^x2+p2 
FD  =  P/2S0 


(10a) 

(10b) 


(10c) 

(lOd) 

(lOe) 
(lOf ) 


'’Piendl,  J.  N.,  "New  Spectral  and  Atomistic  Relations  in  Physics  and  Chemistry 
ol  Solids,"  in  Optical  Properties  of  Solids,  S.  Nudelman  and  S.  S.  Mitra ,  eds. , 
Plenum,  New  York  (1969),  p.  310  ff. 
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The  three  parameters  p,  SQ,  and  x,  appearing  in  the  transition  rate  expression 
are  defined  below,  Equations  (lla-c). 


(11a) 


(lib) 


S 


The  dimensionless  quantity  SQ,  introduced  by  Huang  and  Rhys,®  is  a  measure 
of  the  strength  of  the  anharmonic  coupling  between  a  given  internal  mode  and  the 
acoustic  modes.  As  we  will  see  in  Chapter  4,  S0  is  a  dimensionless  Stokes 
shift,  SQ  =  4fl/2aJD. 

It  remains  to  evaluate  the  internal  mode  reduced  matrix  elements ,  f.- 

#  t  <  ;  Isii 

and/7/K££.  For  this,  we  need  a  specific  model  Hamiltonian  to  be  described  in 
the  next  chapter. 


®Huang,  K.  and  Rhys,  A.,  "Theory  of  Light  Absorption  and  Non-Radiative 
Transitions  in  F-Centers,"  Proc.  Roy.  Soc.,  A204,  406  (1950). 


NSWC  TR  83-90 


CHAPTER  3 

MODEL  HAMILTONIAN  AND  BORN-OPPENHEIMER  SOLUTION 


MOLECULAR  LATTICE  HAMILTONIAN 


We  write  the  Hamiltonian  of  a  molecular  lattice  in  the  following  form: 


H 


n=l 


1 

2 


Mnun  +  V(ux 


• • >u3rN^ 


(12) 


where  un  is  the  n'th  cartesian  coordinate  displacement  from  equilibrium,  N  is 
the  total  number  of  molecules,  r  is  the  number  of  atoms  per  molecule,  and  Mn 
is  the  mass  of  the  atom  associated  with  the  n'th  coordinate.  The  potential 
energy  V  is  a  function  of  all  the  displacements.  If  we  develop  V  in  a  power 
series  we  obtain 


-£*■ 


(2) 

mn 


umun  + 


m,n 


Z(3) 

Vlmn  ulumun  + 

l,m,n, 


(13) 


There  is  no  linear  term  in  Equation  (13)  since  the  un  are  displacements  from 
equilibrium.  Concentrating  on  the  quadratic  term  in  Equation  (13),  we  may  find 
a  transformation  to  real,  dimensionless  normal  coordinates  in  the  form 


3rN 


where  Anj  is  a  real,  orthogonal  matrix.  In  the  new  coordinates,  the 
Hamiltonian  operator  becomes 


(14) 


where 


H  -  S  \  J  (qJ  ■  ^2  V  22  Aijk9i9j ^  +  •• 

J  \  i,j,k 


Z(3)  I  73 

AliAmjAmkvlran  \j —  — - ; — 

1  ,m,n  \MlMmMn^i^i 


i“k 


(15) 


(16) 
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Up  to  this  point,  we  have  made  no  distinction  between  lattice  and  internal 
modes.  Now,  following  Toton,*  we  isolate  one  internal  mode  whose  coordinate  we 
designate  q0.  Greek  subscripts  will  designate  acoustic  mode  coordinates. 

Then,  the  Hamiltonian  may  be  written  in  the  following  form: 


+ 


(17) 


Since  we  are  interested  in  the  interaction  between  the  internal  mode  qQ  and 
the  acoustic  modes  qK ,  we  have  ignored  terms  which  involve  interactions  with 
internal  modes  other  than  q0  and  anharmonic  terms  which  involve  internal  modes 
alone  or  acoustic  modes  alone. 


BORN-OPPENHEIMER  SEPARATION 


Following  the  Born-Oppenheimer  prescription,  we  separate  the  Hamiltonian, 
Equation  (17),  into  the  fast  part 


2  2 
J -  +  (1+2C)  qQ  +  2Bq0  +  2A 

3q2 

o 


(18) 


where 


■ftwoA 


-fiu>0B 

^oc 


fujiKqK 

B<x9Kqx 


and  the  slow  part 


(19) 


*Toton,  E.  T.,  NSWC,  unpublished. 
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where  '►^^(Z)  is  a  normalized  harmonic  oscillator  eigenfunction.  The 
factor  /p  normalizes  the  wave  function  as  a  function  of  qQ.  Both  the  energy 
eigenvalues  and  the  wave  functions  depend  parametrically  on  the  acoustic 
coordinates  q*  through  the  quantities  A,  B,  and  p . 


EQUILIBRIUM  ACOUSTIC  NORMAL  COORDINATES 

The  equilibrium  acoustic  normal  coordinates  are  determined  by  the  energy 
eigenvalues,  Equation  (25),  which  act  as  effective  potentials  for  the  acoustic 
wave  equation.  The  equilibrium  coordinates  are  derived  in  Appendix  F  and  the 
result  [Equation  ( F— 3 ) ]  is 


where  qKS  is  the  equilibrium  value  of  coordinate  q^  when  the  internal 
mode  is  in  state  s. 


(27) 


NON-ADIABATIC  OPERATOR  MATRIX  ELEMENTS 


We  may  substitute  the  wave  functions,  Equation  (26),  into  Equations  (7a, b) 
to  obtain  the  non-adiabatic  operator  reduced  matrix  elements <£K and  £i‘  The 
calculation  is  outlined  in  Appendix  E  and  the  results  are  given  in  Equations 
(E-6)  and  (E-7).  The  matrix  elements  are  functions  of  the  acoustic  coordinates, 
but,  in  the  spirit  of  the  Condon  approximation,  we  will  assume  that  they  are 
constants  with  values  obtained  by  substituting  the  equilibrium  Values  of  the 
acoustic  coordinates  for  the  internal  mode  in  the  initial  state  i.  In  addition, 
we  assume  that  the  third  order  anharmonic  corrections  to  the  Hamiltonian, 
Equation  (17),  are  small  so  that  the  quantities  B,C«1  and  pwl. 


The  results  of  Appendix  E  indicate  that  the  non-radiative  transition  rates 
are  non-zero  only  for  nearest-neighbor  (i+iil)  and  next-nearest-neighbor  (i-*-L±2) 
transitions.  Of  these,  the  nearest  neighbor  transition  rates  will  be  typically 
larger  by  many  orders  of  magnitude.  Thus,  from  Equations  (E-6,  E-7)  we  obtain 


(28) 


(29) 


ESTIMATING  NORMAL  MODE  ANHARMONIC  COUPLING  COEFFICIENTS 

We  are  now  in  a  position  to  determine  the  non-radiative  transition  rate  if 
we  estimate  the  normal  mode  anharmonic  coupling  coefficients  B^  and  CK . 
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Comparing  Equations  (16)  and  (17),  we  can  obtain  expressions  for  CK  and 
in  terms  of  the  third  order  coefficients  of  the  power  series  development 
of  the  potential  energy  of  the  crystal.  In  particular,  we  write  CK  and 
in  the  following  suggestive  forms 


/n 


3k 


M 


3/2 


(30) 


where 


3k 

W72 


^  '  AloAmoAnk  y(3) 
l,m,n,  lmn 


B*X 


V 

l_  V3kX 
/K  m3/2 


jtl L 

'/W^x 


where 


3kX 


/N  M3/2 


^  ^  AloAmtcAnX  y(3) 
l,m,n  /MiMmMn  lmn 


(30a) 


(31) 


(31a) 


In  Equations  (30a)  and  (31a),  M  is  the  mass  of  a  molecule  and  N  is  the  number  of 
molecules.  Since  binding  forces  are  not  long  range,  the  coefficients  V^3^ 
are  significant  only  when  l,m,n  are  nearly  equal.  Therefore,  the  triple  sum  has 
effectively  the  order  of  N  terms.  The  orthogonal  transformation  coefficients 
j  are  typically  sinusoidal  and  proportional  to  N-^2.  Therefore,  CK  is 
proportional  to  N»N-3^2  =  N-^/2.  Apparently,  this  is  true  also  of  but  it  is 

not  clear  that  there  are  not  other  constraints  which  determine  the  non-zero 
values  of  .  It  is  clear,  for  example,  that  the  energy  eigenvalues  in  Equation 
(25)  should  be  essentially  independent  of  the  size  of  the  system  and  so  the 
quantity  B  appearing  in  Equation  (25)  and  defined  after  Equation  (18)  should  be 
essentially  independent  of  N  when  B  is  evaluated  for  q,,-  equal  to  the  equilibrium 
values  qKs.  For  this  to  be  true,  we  require  B^\  to  Be  proportional  to  N-3, 

Then  Equations  (30a)  and  Ola)  define  average  third  order  coefficients  and 

v3kX  * 


Calculations  for  a  linear  diatomic  chain  indicate  that,  for  a  rough 
approximation,  we  may  set  V3K^  ~  f 1  '  *  ( r  e )  and  V  3^.  - —  2“3//2  f''’(re)  where  f(r)  is 
a  pair  potential  describing  the  intermolecular  force  and  re  is  the  equilibrium 
molecular  separation.  Using  this  approximation,  we  may  obtain  an  estimate  of 
the  function  which  appears  in  the  transition  rate  expression  (see  Equation 
(10c)).  The  calculation  is  described  in  Appendix  G.  When  Sc«l  and  Tj)<T<Tp/2S  , 
the  resulting  transition  rate  reduces  to  the  expression  in  Equations  (2,  2a,  2b) 
in  which  we  have  set  Afi  =  ojpSp/2.  The  results  of  Appendix  G  (see  Equation  G-9) 
are  combined  with  the  transition  rate  expressions  (Equations  (10),  (11))  and 
summarized  below  [Equations  ( 3 2— 3 2f ) ] . 
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p 

wi*f  =  ^  WD  —  exp  /g(z0>! — : 

/P  '  '  (flO)  )‘ 


|g(zQ)| 


exp  li¥r  p  +  p"s)  yyp 


(flU)]))  ‘ 


*-[■ 


(i+Dfif^i+l  +  iS 


+  4(2i+l) (P-S)  +  (2i+l)2  P+4 


i6f,i-l  |(2i 


i+1 )2(P~S)2 


.  “fi 


SQ/sinh  (1.  Uficoj)) 

■^P2  +  x2 
SQ  coth  (A.  fifiojj)) 


(32a) 


(32b) 


(32c) 


(32d ) 


(32e) 


(32f ) 


For  up-transitions,  the  final  state  f  =  i+1  so  that  is^  =  m0  where  to0  is  the 
internal  mode  vibrational  frequency.  For  down  transitions,  f  =  i-1  so  that 

“fi  =  -wo- 


ie 
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CHAPTER  4 


RADIATIVE  TRANSITIONS  AND  THE  STOKES  SHIFT 


GOLDEN  RULE  FOR  RADIATIVE  TRANSITIONS 

In  order  to  obtain  the  transition  rates  for  radiative  transitions,  we  can 
adapt  the  formalism  for  non-radiative  transitions  by  replacing  the  non-adiabatic 
operator  <£by  the  dipole  moment  operator  p.  The  expression  for  the  radiative 
transition  rate  corresponding  to  Equation  (4)  is  then 


=  I*  p  ^  ^  Pnj<f  ,m|p|i,n>|  26  (Efm-E£n±fi(D) 


where  p  is  a  factor  involving  the  intensity  of  incident  radiation  (absorption) 
or  density  of  radiation  field  states  (emission)  whose  structure  is  not  important 
for  our  discussion.  The  absorbed  or  emitted  photon  has  frequency  w  with  the 
positive  sign  for  emission  and  the  negative  sign  for  absorption.  If  we  make  use 
of  the  Condon  approximation,  we  can  assume  that  (=<f|pfi>)  is  a  constant  so 

that  the  radiative  transition  rate  can  be  obtained  from  the  last  term  in  the 
non-radiative  transition  rate  expression  by  replacing  with  P  |  W  f i |  ^  (see 

Appendix  C,  Equations  (C-4)(C-9)).  The  result  is  that  the  emission  or 
absorption  at  frequency  <d  is  proportional  to  the  integral 


CD 

(±io )  =  J  e±iwt  +  8(it; 


where  g  is  the  function  defined  in  Appendix  D,  Equation  (D-2). 


HUANG-RHYS  FACTOR  RELATED  TO  MOMENTS  OF  SPECTRAL  DISTRIBUTION 

It  was  first  pointed  out  by  Lax  that  the  Huang-Rhys  factor  is  directly 
related  to  the  moments  of  the  distribution,  Equation  (34)  (see  Reference  4, 
footnote  p.  3).  The  derivation  is  repeated  in  Appendix  H  wherein  it  is  shown 
that  the  Stokes  parameter  Afi  is  the  difference  between  the  means  of  absorption 
and  emission  bands  and  that  Afi  =  2S0<op.  It  is  also  seen  that  the  second 
moment  is  related  to  SQ  (Equation  (H-13)): 


c2  =  up  SQ  coth  (±  B-ficu d ) 
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If  the  distribution  were  gaussian,  it  would  be  possible  to  determine  S0  by 
measuring  the  full  width  at  half  maximum,  Aw,  since,  in  this  case, 

Aw  =  2/2  ln  :  c  (3b) 


SPECTRAL  BAND  SHAPE  ”0T  GAUSSIAN 


Unfortunately,  we  see  by  examining  the  expressions  for  skewness  and  kurtosis 
(Equations  (H-14 ) ( H-15 ) )  that  the  deviation  from  a  gaussian  distribution  is 
considerable  if  S0<<1  which  is  the  case  for  internal  mode  vibrations.  Typical 
internal  mode  frequencies  lie  in  the  range  500  crn-^-  to  2500  cm"^  (10^  sec-^  to 
5  x  10^  sec--*-)  with  typical  infrared  absorption  bandwidths  of  10-50  cm“^. 

Typical  Debye  frequencies  are  around  100  cm-^  (2  x  10 ^  sec-^  or  150°K).  Thus, 
Aw/wq  ~  0.2  which  implies,  from  Equation  (35),  that  S  1/200  and,  therefore,  if 
T>75°K,  that  So<l/200.  This  is  in  sharp  contrast  to  the  situation  for  electronic 
transitions  in  F-centers  where  Aw~2000-5000  cm“^  so  that  So~20  (see  Reference  8, 
footnote,  p.  10). 


HOW  TO  DETERMINE  HUANG-RHYS  FACTOR 

It  is  practically  impossible  to  determine  experimentally  the  moments  of  a 
typical  infrared  spectral  band  due  to  the  presence  of  noise.  Besides  the  mean 
frequency,  the  width  at  half  maximum  is  about  the  only  parameter  which  can  be 
measured  with  any  degree  of  accuracy.  However,  it  is  possible  to  relate  the 
band  width  and  the  Stokes  parameter  even  for  a  non-gaussian  distribution  using 
the  method  described  in  Appendix  I.  A  band  shape  function  with  arbitrary  third 
and  fourth  moments  is  chosen  and  the  moments  are  constrained  to  satisfy 
Equations  (H-13)  -  (H-15).  Then  SQ  can  be  determined  from  the  bandwidth  by 
finding  the  root  of  a  transcendental  equation.  In  the  next  chapter  we  will 
apply  this  procedure  to  the  energetic  material  nitromethane . 


SADDLE  POINT  APPROXIMATION  NOT  GOOD  FOR  RADIATIVE  TRANSITIONS 

By  using  the  method  of  moments,  we  have  not  found  it  necessary  to  actually 
evaluate  the  integral  in  Equation  (34).  This  is  fortunate  since  the  saddle 
point  approximation  which  we  used  for  the  non-radiative  transition  rates  is  not 
a  good  one  for  the  radiative  rates.  This  point  is  discussed  in  Appendix  J. 
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CHAPTER  5 

APPLICATION  TO  NITROMETHANE 


LIQUID  VERSUS  CRYSTALLINE  NITROMETHANE 

In  this  section,  we  apply  the  results  of  the  preceding  sections  to  the 
extensively  studied  and  relatively  simple  condensed  explosive  material, 
nitromethane .  One  caveat  is  necessary,  however.  Under  normal  conditions, 
nitromethane  is  a  liquid  and  even  under  the  extreme  conditions  characteristic  of 
detonations  it  most  likely  retains  the  structure  of  a  liquid.  The  preceding 
results  apply,  strictly  speaking,  to  a  material  which  has  long  range  periodic 
structure.  It  is  not  clear  what  the  absence  of  such  a  structure  would  have  on 
the  predicted  transition  rates.  In  order  to  make  a  direct  comparison  with  the 
results  to  be  presented,  it  will  be  necessary  to  perform  shock  experiments  on 
solid  nitromethane. 


PARAMETERS  TO  BE  DETERMINED 

In  order  to  apply  the  transition  fate  formula  (Equation  (32)),  we  need  four 
parameters  for  a  material.  Two  are  characteristic  of  the  state  of  the  shocked 
material:  the  Debye  frequency  and  the  temperature  T.  Two  are  characteristic 

of  the  internal  mode  whose  transition  rate  is  to  be  determined:  the  vibration 
frequency  <u0  and  the  Huang-Rhys  parameter  SQ. 

The  Debye  frequency  is  a  function  of  compression  which,  in  turn,  is  directly 
related  to  the  shock  pressure  via  the  Hugoniot  relation.  The  internal  mode 
vibrational  frequencies  are  only  slightly  affected  by  compression**  so  that  we 
can  safely  assume  them  to  be  constant. 


HUANG-RHYS  FACTORS  FOR  NITROMETHANE 

The  variation  of  the  Huang-Rhys  factor  SQ  with  compression  is  not  known, 
but  the  bandwidths  of  infrared  absorption  bands  generally  increase  slightly  with 
compression.  In  the  absenbce  of  more  definitive  data  at  the  present  time  we 
will  assume  S0  is  constant.  Table  1  lists  values  of  several  optically  active 
internal  modes  of  nitromethane.  The  values  are  calculated  from  the  bandwidths 


'if 

See  Reference  10,  footnote  p.  23. 

**Typically  less  than  0.02%  per  kilobar  for  the  optically  active  internal 
modes  of  nitromethane. 
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using  che  method  described  in  Appendix  I.  Two  values  ot  SQ  are  given, 
corresponding  Co  Che  cwo  slightly  different  representations  ot  the  band  shape, 
Equations  (1-5)  and  (1-6). 

We  should  note  here  thac  the  optically  active  modes  are  a  subset  ot  cne 
tiny  fraction  of  modes  which  have  zero  wave  number  out  of  the  total  ot 
approximately  Avogadro's  number  of  internal  modes. 


TABLE  1.  HUANG-RHYS  PARAMETER  FOR  SOME  OPTICALLY  ACTIVE  INTERNAL 
MODES  OF  NITROMETHANE  (CH3N02) 


MODE 

iD0(cm”^ ) 

Aw  ( cm- ^  )  * 

SQ (Eq .  1-5) 

S0(Eq.  1-6) 

CH3  rocking  parallel 
to  N02  plane 

1104 

29±6 

1.47xl0-2 

1.39xl0-2 

C-N  stretch 

923 

6 . 5±1 . 1 

l.OlxlO-3 

9.6xl0"4 

N02  symmetric  bending 

663 

19±2. 3 

7 . 3xl0~3 

6.9xl0"3 

NO2  rocking  perp.  to 

N02  plane 

609 

8.211.1 

1 . 58xl0“3 

1.5Uxl0-3 

N02  rocking  parallel 
to  N02  plane 

485 

7.611.5 

1 . 37xl0-3 

1.30xl0~3 

*Bandwidths  were  measured  from  infrared  pbsorption  spectra  obtained  in 
a  diamond-anvil  cell  at  room  temperature  (295°K)  and  pressures  ofi  the 
order  of  5-18  kbar.  Spectra  obtained  by  J.  W.  Brasch,  Jr.  of  NSWC. 


The  data  from  which  SQ  is  calculated  were  obtained  in  a  diamond  anvil 
cell  at  room  temperature  (295°K)  and  relatively  low  pressures  (5  to  18  kbar)  by 
J.  W.  Brasch  of  NSWC.  The  nitromethane  is  in  a  solid  polycrystalline  torm  under 
these  conditions. 

The  values  calculated  for  S0  are  for  a  Debye  frequency  of  2  x  101  rad/sec 
(106.2  cm")  or  152. 7°K).  It  was  determined  that  this  is  a  reasonable  value 
for  the  Debye  frequency  at  standard  temperature  and  pressure  by  taking  the 
centroid  of  the  low  frequency  part  of  the  infrared  absorption  spectrum  of 
nitromethane . 
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TRANSITION  RATES  VERSUS  DEBYE  FREQUENCY  AND  TEMPERATURE 

In  Figure  1,  we  have  plotted  the  transition  rates  given  by  Equation  (32) 
for  the  internal  modes  listed  in  Table  1  using  the  values  of  S0  given  in  the 
next  to  last  column  of  the  Table.  The  transition  rate  is  plotted  as  a  function 
of  u d  for  300°K  (solid  curves)  and  for  2100°K  (dashed  curves).  The  most 
noticeable  feature  of  the  curves  is  the  large  variation  in  transition  rates  with 

relatively  small  changes  in  uQ,  u)q,  and  T.  The  C-N  stretching  mode,  with  a 

frequency  of  1.74x10^  sec-l  has  a  ground  to  first  excited  state  transition 
rate  of  about  10~^  sec-'*'  at  a  temperature  of  300°K  and  top  of  2x10^  sec-'1'  (this 
point  is  beyond  the  scale  of  the  graph).  On  the  other  hand,  the  NO2  symmetric 
bending  mode  with  about  two  thirds  the  C-N  stretching  frequency,  three  times  the 
C-N  stretch  bandwidth,  at  seven  times  higher  temperature  and  two  times  higher 
Debye  frequency  has  a  transition  rate  of  10^  sec-^,  twenty  six  orders  of 
magnitude  greater. 

One  conclusion  to  be  drawn  from  Figure  1  is  that  energy  is  initially  coupled 
most  rapidly  into  the  lowest  frequency  internal  mode,  NO2  rocking  parallel  to 
the  NO2  plane  at  485  cm-*- .  This  is  not  to  say  that  there  are  not  other  internal 

modes  into  which  energy  is  coupled  more  rapidly  since  we  presently  have  no 

transition  rate  information.  We  are  discounting  the  NO2  symmetric  bending 
mode  since  the  calculated  large  transition  rate  is  due  to  its  large  apparent 
bandwidth  and  the  band  may  actually  be  the  superposition  of  two  individual  bands. 

Another  conclusion  to  be  drawn  from  Figure  1  is  that  small  increases  in 
temperature  and  compression  (compression  increases  the  Debye  frequency)  will 
lead  to  large  increases  in  the  rate  of  energy  transfer.  We  also  note  that  the 
relative  importance  of  the  internal  modes  in  energy  relaxation  can  change  with 
changes  in  temperature  and  compression.  For  example,  at  300°K,  the  CH3  rocking 
parallel  to  plane  the  NO2  (CH3  rock  II )  has  a  significantly  lower  excitation 
rate  than  the  NO2  rocking  perpendicular  to  the  NO2  place  (NO2  rock JL ) ,  while  at 
2100°K  the  excitation  rates  are  comparable,  with  the  CH3  rock  II  rate  exceeding  the 
NO2  rock  1_  at  the  higher  values  of  wq. 

Inserting  numerical  values  in  Equation  (2),  we  discover  that  W^  is 
typically  ten  times  greater  than  Wq^,  and  W«3  is  typically  forty  times  greater. 
Thus,  the  relaxation  time  for  energy  distribution  among  the  low  lying  internal 
levels  is  determined  by  the  ground  to  first  excited  state  transition  rate, 

Wqi*  We  cannot  say  anything  about  transitions  between  levels  lying  near  the 
top  of  the  potential  well  since  our  analysis  assumes  a  harmonic  internal  mode 
potential.  At  this  point  we  can  only  say  that,  if  the  transitions  between 
levels  close  to  dissociation  are  also  rapid,  then  the  rate  Wqi  would  be  the 
significant  parameter  determing  the  overall  dissociation  rate.  We  may  expect 
energy  to  be  redistributed  between  internal  modes  more  rapidly  than  it  would  be 
transferred  between  acoustic  and  internal  modes.  Thus,  the  overall  internal 
mode  thermal  relaxation  time  as  well  as  the  overall  dissociation  rate  should  be 
controlled  by  Wqi  for  the  fastest  internal  mode  (which  seems  to  be  the  NO2 
rock  II  in  nitromethane ) . 
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RELATION  OF  DEBYE  FREQUENCY  TO  COMPRESSION 

Since  we  have  repeatedly  stressed  (pun  intended)  the  relation  of  Debye 
frequency  to  compression,  it  will  be  of  interest  to  find  this  relation  for 
nitromethane.  We  can  determine  this  relation  approximately  by  using  the 
expression  for  the  Gruneisen  parameter  which  arises  in  the  Debye  mode  of  a  solid, 


din  u>D 
din  v 


(37) 


where  y  is  the  Gruneisen  parameter  and  v  is  the  specific  volume. ^  Integrating 
Equation  (37)  we  get 


'D  = 


exp 


din  vi 


(38) 


Hardesty  and  Lysne^®  have  calculated  the  thermodynamic  properties  of  shocked 
nitromethane  along  Hugoniot's  for  initial  pressure  of  1  bar  and  initial 
temperatures  of  244°K,  298°K,  and  373°K.  If  we  numerically  integrate  their 
Gruneisen  parameter  data  according  to  Equation  (38),  we  obtain  the  results  shown 
in  Figure  2.  It  is  assumed  that  0)jyo  =  2  x  10^  sec-^.  Figure  3  displays  a  log- 
log  plot  of  the  same  data  as  Figure  2.  From  Figure  3  we  see  that,  above  5  kbar, 
the  curves  are  nearly  linear  indicating  an  approximate  power  law  relation.  The 
dashed  lines  are  isotherms  showing  the  shock  temperature  on  each  Hugoniot. 


TRANSITION  RATES  ALONG  HUGONIOT 

Using  the  result  shown  in  Figures  2  and  3  we  can  plot  the  transition  rates 
versus  pressure  along  a  Hugoniot.  This  is  done  in  Figure  4  for  the  Hugoniot 
with  initial  temperature  298°K.  Again,  Figure  3  is  a  log-log  plot  of  the  same 
data  showing  very  close  to  power  law  curves  above  5  kbar. 


PRESSURE-TIME  CRITICAL  RELATION 

If  we  disregard  the  NO2  symmetric  bending  mode,  the  NO2  rocking  mode 
parallel  to  the  NO2  plane  is  the  mode  whose  transition  rate  is  the  most 
important  in  energy  transfer  from  the  acoustic  modes.  If  the  transition  rates 
calculated  from  Equation  (32)  are  approximately  correct  and  if  the  transition 
rate  for  the  NO2  rocking  (parallel)  mode  is  the  controlling  factor  in  the 
overall  dissociation  chain,  then  the  plot  of  transition  rate  versus  shock 


^Born,  M.  and  Huang,  K. ,  Dynamical  Theory  of  Crystal  Lattices,  (Oxford  Univ. 
Press,  1954),  sec.  4. 

^Hardesty,  D.  R.  and  Lysne,  P.  C.,  Shock  Initiation  and  Detonation  Properties 
of  Homogeneous  Explosives,  Sandia  Laboratories  Report  SLA  74-0165,  May  1974. 
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FIGURE  2.  DEBYE  FREQUENCY  OF  NITROMETHANE  AS  A  FUNCTION  OF  PRESSURE  ALONG  THREE  SHOCK 
HUGONIOrs  WITH  INITIAL  TEMPERATURES  OF  244°K,  298°K,  AND  373°K 
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pressure  for  this  mode  is  a  sort  of  pressure-time  criterion  for  initiation  of 
reactions.  For  it  is  clear  that  a  shock  of  a  given  pressure  must  be  sustained 
for  a  time  which  is  some  multiple  (of  order  unity)  of  Wq[  in  order  for 
significant  dissociation  to  occur. 

COMPARISON  WITH  PRESSURE-TIME  CRITICAL  INITIATION  DATA 

In  this  regard,  de  Longueville,  Fauquignon,  and  Moulard  reported  critical 
initiation  data  in  the  pressure-time  plane  for  several  condensed  explosives 
including  nitromethane.  1  We  have  included  that  data,  shown  as  short  dotted 
line  segments  labelled  DLFM  on  the  graphs  in  Figures  4  and  5.  The  inverse  of 
de  Longueville,  Fauquignon,  and  Moulard' s  time  is  plotted  on  the  ordinate,  Wqi 
It  is  interesting,  though  possibly  coincidental,  that  the  DLFM  data,  over  its 
limited  range,  show  times  that  are  approximately  six  to  seven  times  the  NO2 
rocking  (parallel)  transition  lifetimes  at  the  corresponding  pressures  and  the 
curve  segments  show  roughly  the  same  slope. 


Ude  Longueville,  Y. ,  Fauquignon,  C.,  and  Moulard,  H. ,  Initiation  of  Several 
Condensed  Explosives  by  a  Given  Duration  Shock  Wave,  Sixth  Symposium 
(International)  on  Detonation,  ONR  ACR-221,  Aug  1976,  pp.  105  ff. 
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CHAPTER  6 
CONCLUSION 


We  have  presented  a  quantum  mechanical  calculation  of  the  transition  rates 
for  shock  induced  transitions  between  the  low  lying  internal  molecular  normal 
modes  in  a  molecular  solid.  We  have  assumed  that  the  shock  produces  a 
distribution  of  acoustic  phonons  which  becomes  therraalized  before  any  significant 
internal  mode  phonons  are  created.  This  assumption  seems  to  have  been  justified 
in  the  case  of  nitromethane  in  which  the  shortest  internal  mode  transition 
lifetimes  are -of  the  order  of  nanoseconds  while  lattice  relaxation  times 
determined  from  Van  Vleck's  calculation  are  of  the  order  of  picoseconds  or 
less.  In  particular,  at  the  von  Neumann  spike  pressure  in  nitromethane  (about 
200  kbar),  the  N0£  rocking  (parallel  to  NO2  plane)  mode  has  an  excitation 
time  of  about  4  nsec. 

When  we  compared  the  excitation  lifetimes  with  the  pressure-time  critical 
initiation  data  of  de  Longueville,  Fauquignon,  and  Moulard,  we  found  that  the 
times  were  not  inconsistent  with  the  hypothesis  that  the  overall  dissociation 
rate  limiting  factor  is  the  relaxation  time  (Wq})  for  transferring  energy 
from  the  acoustic  modes  to  a  limited  number  of  interna]  molecular  modes  (NO2 
rocking  parallel  to  NO2  plane  to  nitromethane). 

We  must  reiterate,  however,  that  the  numbers  we  have  obtained  for 
nitromethane  are  subject  to  many  uncertainties,  among  them,  uncertainties  in 
determining  the  Huang-Rhys  factor  Sc  for  each  mode  and  the  uncertainty  in  a 
suitable  choice  of  Debye  frequency  10 p.  Small  changes  in  both  of  these 
quantities  lead  to  large  changes  in  the  transition  rates.  Also,  for  the  great 
majority  of  modes  which  are  not  optically  active,  it  is  not  possible  to 
determine  SQ.  Perhaps  the  results  of  neutron  scattering  experiments  may  give 
useful  information  on  these  modes.  Other  problems  are  the  question  of  the 
validity  of  the  Condon  approximation  [Equation  (9)],  the  determination  of  the 
anharmonic  coupling  coefficents  [Equation  (30)  and  (31)],  the  validity  of 
isolating  one  internal  mode,  neglecting  the  interaction  between  internal  modes, 
and  the  use  of  a  single  frequency  wp  to  characterize  the  acoustic  spectrum. 
Finally,  the  shock  data  available  for  nitromethane  is  for  the  liquid  state, 
whereas,  the  calculations,  strictly  speaking,  apply  to  the  solid  state.  We  hope 
that  experimenters  will  be  encouraged  to  undertake  shock  experiments  on  solid 
nitromethane  in  order  to  obtain  both  Hugoniot  data  and  critical  initiation 
data.  In  order  to  apply  the  results  presented  here  to  other  solid  explosives  it 
is  necessary  that  the  explosive  be  homogeneous  and  have  an  internal  mode 
spectrum  clearly  distinguished  and  well  separated  from  the  lattice  spectrum. 
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APPENDIX  A 


BORN-OPPENHEIMER  APPROXIMATION 


Suppose  we  can  split  the  coordinates  of  a  system  into  two  groups,  a  fast 
subsystem  {r£}  and  a  slow  subsystem  {q<}*  The  Hamiltonian  of  the  system  is 

HCr^q,,.)  =  Tr  +  Tq  +  VCr^.q*)  (A-l 

where  Tr  and  Tq  are  the  kinetic  energy  operators  for  the  fast  and  slow 
coordinates,  respectively.  Let 

Hr(qK)  =  Tr  +  V(ri>q<)  (A~2 

and  consider  the  solutions  i(is  of  the  eigenvalue  problem 

Hr^s^ij^ic  )  =  e s ^ 9ic  ^  ^s(ri»q<  )  (A-3 

The  eigenvalues  es  and  the  wave  func t ions.  i()  s (r £ )  depend  parametrically  on  the 
coordinates  q^. .  Suppose  the  set  { ip  s C r  ^ , q^. ) >  is  a  complete  orthonormal  set  in 
r-space  for  every  of  values  of  q^. .  Then  any  solution  ¥  of  the  original  system 
Schrodinger  equation, 

H¥  =  E¥  (A-4 

can  be  written 

^i’^  =  ♦s  (ri>9K)  (A-5 

s 

Then  Equation  (A-4)  is  equivalent  to 

^[es(qK)  +  Tq  -  Ej^s(q<)  *s  (r^q*)  -  0  (A-6 

s 

Now,  consider  Tq  and  'l's^ri>9<^  as  operators  in  q-space.  Then 

Tq$  s  ^  ri  >  )  “  t  "^q  j't' s  ^ri  >  ^  J  +  't,s^ri,9ic  ^  Tq  (A— 7 

Substituting  (.A-l)  into  (A-6)  we  obtain  the  following  equivalent  version  of 
(A-4): 

(}£+<£)  ¥  =  E¥  (A-8 
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where 


=  X]  4's(ri’c’K)  Tq  +  es(<«K)[ 
s 


(A-9a) 


(A-9b) 


In  Equations  (A-9a,b),  is  the  adiabatic  Hamiltonian  and  <£  is  the  non-adiabatic 
operator.  The  usual  Born-Oppenheimer  basis  is  obtained  by  ignoring  oC  and  solving 
the  equation  V  =  E'i ,  or,  equivalently, 


|  ^q  +  es^(l<^|  ^sn^ie^  ^sn^sn^  ^ 


(A-10) 


where  s,n  are  the  fast  and  slow  subsystem  quantum  numbers,  respectively.  The 
Born-Oppenheimer  basis  states  are 

*sn^ri *9*  ^  =  ^s^ri»^<^sn^9ie) 


( A- 1 1 ) 
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APPENDIX  B 

NON-AD IABATIC  OPERATOR 


Equation  (A-9b)  defines  the  non-adiabatic  operator.  What  does  this  operator 
look  like  if  the  slow  subsystem  consists  of  a  set  of  harmonic  oscillators?  In 
this  case  the  kinetic  energy  operator  Tq  is 

Tq  V 

tc 

In  Equation  (B-l),  the  kinetic  energy  has  been  written  using  dimensionless 
coordinates  q  =  /mu /ft  x  and  the  normal  mode  frequencies  are  denoted  by  coK.  The 
action  of  the  non-adiabatic  operator  on  fsn  is  therefore 


2 


(B-l) 
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APPENDIX  C 

EVALUATION  OF  SUM  OVER  STATES 


Introducing  the  integral  representation  of  the  delta  function,  Equation  (4) 
becomes 


CD  _ 

W(i-*-f)  =  i_  /  /  exp 

*2  A  ki 


2 

|i(Efm-Ein)t/-fi|  Pn|<fm|<£|in>|  dt  (C-l) 


In  the  lowest  order,  the  acoustic  mode  wave  functions  are  products  of  harmonic 
oscillator  wave  functions: 


►in  -n 


U) 

►in 


(C-2) 


The  energy  levels  are 


Ein  =  Ei  + 


?M) 


-flcu  „ 


and  the  square  of  the  matrix  element  appearing  in  (C-l)  is 


,  .  (k) 

2 

|<fm|£|in>|2  =  yjlK{i\ 

K 

(<)  i3Hn 

4fn  a  *> 

K  1  3qK 

n 

v/< 

(p)  |  (p) 

^ fm  ♦in  > 

p  1  P 

(C-3) 


U) 


V"  J  f*  (O.^in  (X)  3(>,ink 

+  tcfi*Scfi  <4>fr  1  -w,kr-  I  x 


ic  ,  X 

ic/x 


3c5k  X 


(x) 

FT  (v) 

>11  <*fm 
3(3X  p/<  v 


(p)  TT  (v)  (v) 

♦  in  >  1  -*-<(♦ fm  '►in  > 

v/X  v'  v 


*  (£' O  YAn  i  <£’  i .£V  n 

\X  ’  <  tc  3qK  k  k  U=|C 


(p)  .  (p) 

<►  fm  |  4>in  > 


+  c.c. 


C-l 
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y Td  2  tt  .  <v> 

+  <fi  ■*■  Ofm  ^in  ^ 


IC-4; 


M  IJ 


Substituting  (C-4)  in  (C-l)  we  obtain 


CD 

W(i+f)  =  i-  (  exp(iwfit)  ('Ffi(t)  +'Efi(t)  +HfL(t)  +Mfi(t)l  at 

-s  .  1  '  (c- 


where 


Wfi  -  (Ef-E^) //fi 


(C-5) 


(C-5a) 


Ffi(c)  FI  Gpfi(t' 


(C-5b) 


Efi^^  *  y  '°4fiAfi  ^Kfi^c^Exfi^-t^  rr  cyf^t) 


(C-5c) 


^fi^  “  VAfi  n  Gufj(t )  y  SJ^-\ f  i  +  complex  conjugate 

x  y^x  X  (C-5d) 


Mfi^c^  ^  j^Kii  I  lGpf  i  ^ c - 


(C-5e) 


.  ,  i(m-n)w  t  j  (y)  (y)  I  2 

■  -  /  4  y  j  Pn  e  pC^fm  |  ♦  in 


m=o  n=o 


(C-6a) 


00  00 

Fxfi(t)  =  ^  ^  pn 


m=o  n=o 


i(m-n)<D  t  I  (x  )  „  (x)  2 

,  K  <4»fra  |-L_*in> 

1  3qk 


(C-6b ) 
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^  tt— \  \  i(m-n)u)  t  (< )  (<  )  (<)  (<) 

f  i ( c )  ~  /  v  e  <  fm  |  ~ —  $  in  '  fm  I  ♦  in  )■  ( C-6c ) 

m=o  n=o  .  4|C 

The  sums  in  Equations  (C-6a,b,c)  can  be  evaluated  using  the  expression  for  the 
density  matrix  of  the  harmonic  oscillator  in  the  coordinate  representation 
(Slater's  sum  or  Mehler's  formula). The  results  are 


^icfi^^  ~  exP  {  ^^icf  i| 

(C-7a) 

Ftcfi^c^  ~  ^  +  [^Kfi] 

(C-7b) 

f  i ^ ^  =  E|C(t)^icfi  GK(t) 

(C-7c) 

where 

cosh(io)  t  +  T  Rfiw  ) 

FK  ( t )  =  <  2  < 

2  sinh(A.  Bfiuj,,.  ) 

(C-8a) 

EK(t)  =  F^Ct)  -A.  coth  (A.  Bfiu)K ) 

(C-8b) 

Aicfi  =  qxf  “  ^xi 

(C-8 c) 

In  Equation  (C-8c),  qKf  and  qK£  are  the  equilibrium  values  of  the  acoustic 
coordinates  when  the  fast  subsystem  is  in  states  f  and  i,  respectively.  Using 
the  results  above,  we  can  write  the  expression  for  the  transition  rate  as 
follows : 

CO 

W(i-*-f)  =  A_  (  exp(itufit)Ff j.(t)Gf£(t)dt 
fi2  'i 

(C-9) 

where 

~  ^  y  1  f  i  1  ^K(t)  +  ^  f  j  ^<fi  E,(t)^  Xfi  AXfi 

K  K  X 

Ex(t) 

+  Wfi  ^°4f  i  ^xfi  +  Mfi  ^  Sj^-k  fi  ^xfi 

X  X 

^“^See , 

for  example,  O'Rourke,  R.  C.,  Phvs .  Rev.,  91,  263  (1953). 

C-3 
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+ 


‘fi 


I>Kfi 


(C-9a) 


(C-9b ) 


(C-9c) 


S 


—  coth 
2 


(i.  lMiulc  )Ak  f  £ 


(C-9d) 


C-4 
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APPENDIX  D 


SADDLE  POINT  INTEGRAL  FOR  TRANSITION  RATE 


It  is  convenient  to  make  the  change  of  variable  z  *  it  in  the  expression 
for  W(i-*f)  in  Equation  (9).  Then  Equation  (10)  can  be  written  in  the  form 


W  ( i-*-  f  ) 


=  — L-  f  Ffi(-iz)  e§^z^dz 


(D-l ) 


where 


Zcosh(u  z+-r  JB'fuo  )  2 

- - —  A*  fi  -  S 

K  2  sinhCi.  ) 


(D-2) 


At  a  saddle  point,  g'(z0)  *  0  so  that 


“fi  +  ^  “ 


sinh(w  z+-r  Bfiu)  )  2 

x  2  x  A 

x  - ; - '  ~  '  °xfi 

2  sinh(i  13fvwK  ) 


(D-3) 


Equation  (D-3)  has  many  complex  roots  but  precisely  one  real  root.  There  is  at 
least  one  real  root  since  the  hyperbolic  sine  function  goes  to  +»  when  z*±». 
There  cannot  be  more  than  one  real  root  because  g"(z)  is  never  zero  for  real  z 
implying  that  g(z)  is  monotonic  [See  Equation  (D-4)]. 


g"(z) 


2  cosh(w  z+-r-(3fiu)  )  2 

0K - £ — l - £_  tKi 

2  sinh(i.  (SfiuiK) 


(D-4) 


In  particular,  g"(zo)>0  so  the  direction  of  steepest  descent  is  along  the 
imaginary  axis.  Near  the  saddle  point 

g(z)  a  g ( z o )  +  i  (z-z0)2g"(z0) 


(D-5) 
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Thus,  deforming  the  contour  of  integration  to  pass  through  the  saddle  point,  we 
get 


00 

i-*-f )  =  —  Ffi(-iz0)exp  (g(z0)}  f  exp  {-  -  t2  g"(z0)}dt 

J  2 


and  the  final  result  is 


W(i-*-f)  =  —  -*/ — —  Ff£(-iz0)  exp  { g C z0 )} 

*2  *g"(z0) 


(D-6) 


(D-7) 


More  details  on  the  saddle  point  integral  may  be  found  in  Perlin's  review  article 
(Reference  3,  footnote,  p.  3). 

If  we  can  replace  the  acoustic  frequencies  by  a  single  average  frequency 
up,  a  considerable  simplification  is  achieved.  Then 

Ffi<-izo)  |^efi  |  FD*|Z^fi  |FD  *  \  coth^y 

+  ]p4fiAKfi  +^^Xfi  ‘  \  COth^  ^ 


Y  ^<fi 


(D-8) 


where 


cosh(u-  z  +  (Jfiuj-.) 
_ Do  z _ D_ 

2  sinh(i  frfiwj)) 


Equation  (D-3)  which  determines  z0  becomes 
P  =  x  sinh  (tun  z^  +  i  ftfiwn) 


(D-8a) 


(D-91 
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APPENDIX  E 

EVALUATION  OF  ^fi  AND  W  fi 


Using  the  result  in  Equation  (26)  we  obtain 


(*£'  §)’/ 


(HO) 

/p"  f  (Z) 


—  I 

3<1<  ( 


(HO)  f 
/p"  H  (Z)  >  dqQ 


(E-l) 


and,  transforming  the  variable  of  integration  from  qQ  to  Z,  the  result  is 

(HO) 

/  3*.\  f  (HO)  d*.  (Z)  /  \ 

(*f,  -^)  =  /  *f  (2)  - -  f5-  dZ  +  -i-  (in  /p  )  6fi 

V  3q</  •/  dZ  3qK  3qK  \  / 


(E-2) 


Since  f*i,  the  second  term  in  Equation  (E-2)  is  zero.  Differentiating  Equation 
(23)  we  obtain 


flU), 


3Z_ 
8  9k 


\ p-4  c<z  +  2p-3  ***<&  ‘  2p"7  B  ^ 


(E-3) 


Using  Equations  (E-3),  (E-2),  and  (7a)  wt  obtain 


4fi  =  -^-Up-4  c*<flzf7l 1 


fw0  )2 


i> 


+  2 


“  P“7  BC^  <f|i^|  ii 


(E-4) 


where  <f'Z  3/3Z  i>  and  <f|3/3Zji>  are  the  usual  harmonic  oscillator  matrix 
elements.  These  may  be  evaluated  most  conveniently  by  writing  Z  and  3/3Z  in 
terms  of  the  creation  and  annihilation  operators,  a+  and  a,  and  then  using  the 
raising  and  lowering  properties  of  these  operators. 


Z  -  I  (a  +  a+) 

fl 

i-  =  I  (a  -  a+) 

3Z  /2 


(E-5a) 

( E-5b ) 


E-l 
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a>  =  /n+T  n  + 1  > 


n>  =  /n  n~l> 


(E-3c ) 
(E-5d ) 


The  end  result  is 


*4cfi  =  ~  |^P"3  BicX9X  ~  P“7  bckJ  [/T  <S  f ,  i  —  1  *  /i+1  6f  ,i  +  lj 

+  I  p-4  CK  pi(i-l)  6  f  ,1-2  ~6  f  £  -  /(i+1 )  (i+2)  6tji+2jj  (E-f 


A  similar,  though  lengthier,  calculation  yields  Since  the  result  has 

many  terms,  we  give  here  only  the  lowest  order  term  involving  powers  and 
products  of  C<  and  B^x : 


^ic  f  i  2  ^  ^  f ,  i  - 


1  -  /i+l  6f 


:,  i  +  1^ 


(E-7) 


E-2 
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APPENDIX  F 

ACOUSTIC  MODE  WAVE  EQUATION  APPROXIMATE  SOLUTIONS 


The  internal  mode  eigenvalues,  Equation  (25),  are  the  effective  potentials 
which  determine  the  acoustic  mode  Born-Oppenheimer  wave  functions  via  Equation 
(22).  To  obtain  an  approximate  solution,  expand  the  eigenvalues  e f  in  powers  of 
q<  up  to  quadratic  terms.  This  is  a  good  approximation  if  the  anharmonic 
corrections  to  the  total  Hamiltonian  are  small  so  that  B,C<<1.  Then, 


(s+i)  CK 
2 


+  -f^q* 


(F-l ) 


and 


32c  .  C  C 

_ k 6 kX  *(s+i)  JL2. 

3qK3qK  2  f*>0  ' 


Thus,  the  equilibrium  acoustic  coordinates  are 


qKS  *  -(s+I)  -~— 

2 


(F-2) 


(F-3) 


and  the  effective  potential  is 

€s  ~  £s  ^s  ^  +  s  ^  ^ 

K 

In  Equation  (F-A)  we  have  ignored  the  off-diagonal  quadratic  terms  since  they  are 
small  under  our  previous  assumptions.  These  terms  lead  to  mixing  of  the  acoustic 
coordinates  and  subsequent  modifications  of  the  acoustic  frequencies.  In  terms 
of  the  variables  q^  =  qK  -  qKS,  the  acoustic  wave  equation,  Equation  (22), 
becomes 


* 


'  sn 


=  Lsn$ 


sn^sn 


(F-5) 


F-l 
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where 


J s  =  (s+I)-fio)0  -  1  -ftuiic  q<  s 


( F-6 ) 


Thus  the  total  energy  eigenvalues  are 

Y  inK+^^K 


Esn  + 


(F-7) 


F-2 
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APPENDIX  G 

APPROXIMATE  TRANSITION  RATE  FORMULA 


We  approximate  the  coefficients  and  B >  using  Equations  (30)  and  (31)  as 
our  guide: 


fxB 


/N 


UD  w0 


3/2 


f2fl 


N  uD 


3/2 


(G-l) 

(G-2) 


We  have  denoted  1/2  f"(re)  by  B.  The  quantities  f^  and  f2  are  dimensionless 
numbers  which  we  will  later  set  to  \//T  and  2,  respectively.  Substitute 
Equation  (27)  into  (9e)  and  substitute  the  result  into  Equation  (lib)  to  obtain 
the  Huang-Rhys  factor  for  a  transition  from  state  i  to  state  f: 


So 


ifcj)2  £ 

2-fi2w2  V 
D 


2 


(G-3) 


The  sum  contains  N  terms,  so  upon  substituting  Equation  (G-l)  into  (G-3)  we  get 


So  = 


1 

2 


(G-4) 


In  the  following,  we  set  f=i+l.  Using  Equations  (28),  (9e),  (27),  (G-l),  (G-2), 
and  (G-4)  we  obtain 


K,i+l,i 


4|_2.)  (i+l)(2i  +  l)2  soWiu>j))2 

V 1/  “o 


jl  +  (2i+l)S0w 


(G-5) 


G-l 


.  j 

.  I 
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The  second  term  in  curly  brackets  in  Equation  (G-5)  arises  from  the  second  term 
in  Equation  (28).  In  our  application,  SQ  is  much  less  than  one,  therefore,  we 
can  omit  the  second  term  in  Equation  (28).  Similarly,  we  obtain 


i*l>il  ’2fe)  (i+l)(2i+l)2  S0(-hto  D)2  u D/o)0 


(G-6) 


E^-,i*i,i  YA.i+i-1  Ax  =  8(-r\  (i+1)(2i+1)  s0  (^d)2  u)D/u>0 
V  x  vv  (g-7) 


I  y]  ^tc ,  i+1 ,  i  I  =  (i+1)  s°  “D^o 


(G-8 ) 


Substituting  Equation  (G-5)-(G-8)  into  Equation  (D-8)  we  obtain 

(fV 


Fi+i ,  i(-izQ)  =/_i j  (i+l)(fkoD)2  Sq  |(2i+l)2  [P+(P-S)2] 


+  4(2i+l ) (P-S)  +  4 


} 


(G-9) 


where  P  and  S  are  defined  in  Equations  (lOd)  and  (lOe).  If  we  now  substitute 
Equation  (G-9)  into  Equation  (10),  we  obtain  an  expression  for  W(i-*i+l) 
depending  on  T,  u>0,  u>D,  and  SQ.  All  the  third  order  anharmonic  coupling  effects 
are  condensed  into  the  Huang-Rhys  factor  SQ. 


If  SQ«1  and  the  temperature  is  not  too  large,  we  may  simplify  the  result 
further.  The  upper  temperature  limitation  is  determined  by  the  condition  x2«l 
where  x  is  defined  in  Equation  (11c).  Thus  the  condition  is 


L  I3fi wD  >  arcsinh  SQ  =  SQ 


(G-10) 


Rearranging  (G-10)  we  obtain  the  condition 
T  <  T0/2So 


(G-ll) 


where  Tq  is  the  Debye  temperature.  Under  this  condition,  S« 1  and  P»  |p|  .  If  we 
set  f2  =2,  =  1//2",  and  Afl  =  u>D  SQ/2,  we  obtain  the  result  given  in  Equations 
(2,  2a,  2b). 


G-2 


I 


'1 

■’4 


1 


I 


* « 
- 


*.  i 

[4 


■U 


.3 

: 

•4 


§ 


■  i 

ft 

1 
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APPENDIX  H 


MOMENTS  OF  SPECTRAL  BAND  SHAPE  FUNCTION 


The  r'th  moment  of  the  distribution  I  (u> )  is  defined  by 


<u>r>  = 


—00 

/oo 

I (u)  )du> 


(H-l) 


where  I(w)  is  defined  in  Equation  (34).  In  order  to  evaluate  <ojr>  note 
that 


re  t  dw  =  2ir ir  d  6  (t ) 
dt1' 


where  6(t)  is  the  Dirac  delta  function.  Thus 

J l(u))o)rdu)  =  2ir  ir  d.  ^c2  dt 


After  r  integrations  by  parts  we  obtain 


I  (oo  )(D 


rdw  =  2ir  .  - -d.!_  r.«<it)l 

d(it)r  L  J  t=c 


From  Equation  (D-2)  we  note  that  g(0)  =  0  so  that 


<wr>  =  (±I)r  .dI..el(.2) 


(H-2) 


(H-3) 


(H-4) 


(H-5) 


The  positive  sign  is  for  absorption  and  the  negative  sign  is  for  emission. 
The  average  absorption  or  emission  frequency  is 


w  =  <u)>  =  tg 1  (0 )  - 


;  (w  f  i  +  2.  Aft  ) 


(H-to) 


H-l 
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where 


a« 


■E 


(H-7) 


For  absorption,  let  w0  =  <D££>0  and  for  emission  let  u)0  =  -u)££>0.  Then 


a)  =  u)n  +  ±  Afi 
°  ~  2 


(H-8) 


Thus,  the  Stokes  parameter  Afi  is  the  difference  between  the  mean  absorption  and 
emission  frequency.  The  second  moment  of.  the  distribution  is 


j2=<(u)-io)2>  =  g"(0)  =  L  ^(ii^  coth  (i  &RtDK)AK 


and  the  third  and  fourth  moments  are 


<(ii)~(i))^>-±.g,,,(0)  ^  ^  0),,. 


3  2 
A- 


(H-9) 


(H-10) 


<(a)-m)4>  =■  g<4>(0)  +  3[g"(0)]2 


EA  2 

(oK  coth (i  +  3[g"(0)]2 

K 


(H-ll) 


Again,  introducing  an  average  acoustic  frequency  wj),  we  can  express  the 
moments  in  terms  of  the  Huang-Rhys  factor: 


Afi  -  2S0toj) 
2 


2  _ 

a*  -  too 


SQ  coth  (A  frfujn) 
2 


Y3  =  coth  (1.  Sfiaij) 

Y4  =  <(qd±>l>  -  3  =  [so  coth  (I  efm>D) 

L  2  J 


-1 


(H-12) 

(H-13) 

(H-14) 

(H-15) 


The  skewness  Y3  and  kurtosis  Y4  are  measures  of  the  deviation  from  a  gaussian 
distribution . 


H-2 
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APPENDIX  I 


DETERMINATION  OF  STOKES  PARAMETER  FROM  BANDWIDTH 


A  straightforward  method  of  obtaining  a  distribution  with  a  finite  number 
of  specified  moments  would  be  to  take  the  product  of  a  polynomial  with  the  error 
function,  choosing  the  coefficients  of  the  polynomial  appropriately.  An 
equivalent,  but  easier,  approach  is  to  take  a  linear  combination  of  the  error 
function  and  its  derivatives  with  appropriate  coefficients: 


I(x)  =  an4>(n)(x) 

n=0 


where 


4>  (x)  =  exp 


|-HA 


(i-D 


(1-2) 


Integrating  by  parts  we  obtain  the  formulas 


f  x2r  l(x)dx  =  f- ?J.) 
•L  or 


r  „n 
r— '  2  a 


2n 

TF^TT 


(i-3) 


f  x2r+1  I (x)dx  =  Y  l_!2n±; 

J  or  Tr-nTT 


(1-4) 


If  we  wish  I(x)  to  be  a  normalized  distribution  with  zero  mean,  unit  standard 
deviation,  skewness  Y3,  and  excess  y 4,  then,  from  Equations  (1-3)  and  (1-4)  we 
find  that  aQ  =  1,  a^  =  0,  a2  =0,  a^  ~  -73/6,  and  a^  =  -y^/24.  Thus 


I(x)  =  *(x)  -  13  *<3 >(x)  +  I*  *(4>(x) 

3  J  4! 


(1-5) 


1-1 
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Now  let  x  =  (w-u))/a  and  substitute  for  o,  Y3>  and  Y4  from  Equations  (H-13), 
(H-14),  and  ( H— 15).  If  we  then  set  I(x)  =  I(0)/2,  we  obtain  an  equation  which 
can  be  solved  numerically  for  S0,  given  (io-uO/io j).  Another  distribution,  quoted 
by  Lax  (Reference  5,  footnote  p.  4)  as  being  due  to  Edgeworth,  includes  a  sixth 
derivat ive 

I(x)  -♦(x)  =  ^3)(X)  +  14  ^  C4>  (x)  +  4)  (6)(x)  (1-6) 

The  distribution  (1-6)  provides  a  not  significantly  better  match  for  the  higher 
moments . 


1-2 
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APPENDIX  J 

VALIDITY  OF  THE  SADDLE  POINT  APPROXIMATION 


The  saddle  point  approximation  for  the  transition  rate,  outlined  in  Appendix 
D,  is  the  first  term  in  an  asymptotic  expansion  in  the  sense  defined  by 
Poincare.3  In  order  for  the  first  term  to  be  a  good  approximation  to  the 
integral,  the  second  term  must  be  much  smaller.  The  expansion  for  an  integral 
of  the  form  in  Equation  (D-l)  has  been  worked  out  by  Hoare  for  the  case  when  the 
function  g(z)  and  its  derivatives  are  of  the  order  of  N  where  N  is  the  asymptotic 
expansion  parameter  and  is  presumed  large. In  our  case,  the  relevant 
asymptotic  parameter  is  p,  defined  in  Equation  (D-9 a)  in  the  average  acoustic 
frequency  approximation.  The  ratio  of  the  second  term  in  the  asymptotic 
expansion  to  the  first  term  must  be  much  less  than  one.  Using  Hoare 's  result, 
we  obtain 


l 

8 


(4) ,  . 

S  (2o) 

[g"(z0)]2 


5[g(3)(zo)]2 

3(g"(z0j]3 


(J-1) 


In  the  average  acoustic  frequency  approximation,  the  derivatives  of  g  are: 


g(2n) (zQ) 


2n 

w  P 
D 


n  >  1 


(J-2) 


where  P 


g( 2n+l ) (zQ)  =  UD 


2n+l 


;  n  ^  1 

(J-2)  and  (J-3)  in  (J-1),  we  get 


(J-3) 


(J-4) 


J-ljef freys ,  H.  and  Jeffreys,  B.  S.,  Methods  of  Mathematical  Physics,  3rd  ed., 
Cambridge  University  Press,  1956,  c.  17. 

J-2noare>  m.  R. ,  J.  Chem.  Phys.  52 ,  5695  (1970). 
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From  (J-4)  we  see  that  the  second  term  in  the  expansion  is  less  than  1/12  the 
first  term  if  |p|  is  greater  than  one,  independently  of  the  value  of  x 
[=  SG/sinh  (3fia)p/2) .  We  also  see  that  the  second  term  is  less  than  1/8  the 
first  terra  if  x  is  greater  than  one.  Tn  fact  the  expansion  is  asymptotic  in  the 
parameter  P  so  that  the  approximation  becomes  better  if  either  x  or  p  becomes 
large . 

In  the  case  of  non-radiative  transitions,  | pf > 1 ,  so  the  criterion  is 
satisfied  whatever  the  value  of  SQ.  In  fact  Jp|»l  is  required  for  the  Born- 
Oppenheimer  approximation  to  be  valid. 

For  radiative  transitions,  however,  the  saddle  point  approximation  is  not 
good.  For  radiative  transitions, 

p  =  -  (<»f£  ±  w)/wD  (j-5) 

where  the  positive  sign  is  for  emission  and  the  negative  sign  is  for  absorption 
at  frequency  w .  The  range  of  values  of  p  which  are  of  interest  (that  is, 
values  corresponding  to  frequencies  within  the  spectral  band)  is  then 


where  Aw  is  the  bandwidth.  Since  SQ  typically  is  much  less  than  one,  both 
|p|  and  x  will  be  much  less  than  one  so  that  condition  (J-4)  cannot  be  satisfied 
for  values  of  p  which  correspond  to  frequencies  lying  within  a  spectral  band. 
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